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ABSTRACT 

We compare the predictions of three physical models for the origin of the hot halo gas with the 
observed halo X-ray emission, derived from 26 high-latitude XMM-Newton observations of the soft 
X-ray background between I = 120° and I = 240°. These observations were chosen from a much 
larger set of observations as they are expected to be the least contaminated by solar wind charge 
exchange emission. We characterize the halo emission in the XMM-Newton band with a single- 
temperature plasma model. We find that the observed halo temperature is fairly constant across the 
sky (~(1.8-2.4) x 10^ K), whereas the halo emission measure varies by an order of magnitude (^^0.0005- 
0.006 cm~^ pc). When we compare our observations with the model predictions, we find that most of 
the hot gas observed with XMM-Newton does not reside in isolated extraplanar supernova remnants 
- this model predicts emission an order of magnitude too faint. A model of a supernova-driven 
interstellar medium, including the flow of hot gas from the disk into the halo in a galactic fountain, 
gives good agreement with the observed 0.4-2.0 keV surface brightness. This model overpredicts 
the halo X-ray temperature by a factor of ~2, but there are a several possible explanations for this 
discrepancy. We therefore conclude that a major (possibly dominant) contributor to the halo X-ray 
emission observed with XMM-Newton is a fountain of hot gas driven into the halo by disk supernovae. 
However, we cannot rule out the possibility that the extended hot halo of accreted material predicted 
by disk galaxy formation models also contributes to the emission. 
Subject headings: Galaxy: halo — ISM: structure — X-rays: diffuse background — X-rays: ISM 



1. INTRODUCTION 

Observations of the diffuse soft X-ray background 
(SXRB) indicate the presence of ~(l-3) x 10^ K X-ray- 
emitting gas in the interstellar medium (ISM) of our 
Galaxy. Early observations with rocket-borne instru- 
ments led to the conclusion that the diffuse 1/4- keV 
emission was dominated by emission from ^^1 x 10^ K 
plasma in the Local Bubble (LB), a ~10 0-pc cavity in 
the ISM in which the Solar System resides (| Sanders et al.l 
Il977t iSnowden et al.lll990[ ). The discovery of shadows in 
the 1/4-keV background with ROSAT showed that there 
was also gas with T '--^ 1 x 10 ^ K beyond the Galac- 
tic disk, in the Galactic halo (jBurrows fc Mendenhalll 
llQQll : ISnowden et aLllTOQll ) . Higher-energy emission data 
from ROSAT, XMM-Newton, and Suzaku, and X-ray ab- 
sorption line data from Chandra, show the presence of 
hotter gas in the Galactic halo, with te mperatures up 
to ^ 3 X 10*^ K (Kuntz & Snowdcn 200(1; lYao fc Wan 
20051 [2007 : Smith ct al. 2007; Galeazzi et al. 200 



Henlev fc Shelton.,2008; ,Yao et al..,2009; ,Lei et al... 20091 : 
Yoshino et al.ll2009D T 



Several possible sources for the hot halo gas have 
been suggested, including supernova- (SN) and stel- 
lar wind-driven outflows from the Galactic disk (e.g. 



Shapiro fc FieIdllT976l:lB^egmanllT980HNorman fc Ikeuchil 
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|dbh@ physast.uga.edu| 

'TDepartment of Physics and Astronomy, University of Geor- 
gia, Athens, GA 30602 

^ Department of Astronomy, Columbia University, 550 West 
120th Street, New York, NY 10027 

^ Department of Astrophysics, American Museum of Natural 
History, 79th Street at Central Park West, New York, NY 10024 



material (predicted by s i mulations of disk galax y for- 
mation; iToft et al.l 120021: iRasmussen et a l.' '2009'), and 
in situ heating by ex traplanar SNe ([Shelton 200(| 
iHenlev fc SheltonI l2009f ) . X-ray spectroscopy is essen- 
tial for determining which process or processes have pro- 
duced the ~(l-3) X 10® K gas in the halo. In princi- 
ple, the observed ionization state could be used to dis- 
tinguish between the different models. For example, 
gas heated by SNe could be underionized if heated re- 
cently, or ove rionized if heated in the distant past (e.g., 
lSheltonlll999f ). and gas that has recently burst out of 
the disk, cooling rapidly, will be dra stically overion- 
ized (jBreitschwerdt fc Schmutziei1ll994[ ). The elemental 
abundance ratios could also, in principle, be used to dis- 
tinguish between models, as the abundance pattern of 
the hot gas may depend on whether it is of Galactic or 
extragalactic origin. 

In practice, it is not easy to use arguments based 
on the ionization state or the abundances to distin- 
guish between models, as collisional ionization equi- 
librium (CIE) models with solar abundances generally 
provide good fits to the observed X-ray sp e ctra (e.g. 
Galeazzi et al. 2001^ IHenlev fc SheltonI 120081 : iLei et al.l 
2009; Yosh ino et all l2009(). although supersolar [Ne/0] 
and [Fe/0] abundance ratios have been reported for some 
sightlines (Yo shino etHI 120091: lYao et al.l |2009[) . Here, 
we use a different approach. We fit CIE models to 26 
XMM-Newton spectra of the SXRB, obtained from ob- 
servations between I = 120° and I = 240° and with 
|6| > 30°. These fits yield temperatures and emission 
measures for the halo. We then compare the measured 
distributions of these quantities to those predicted by 
two physical models of the hot halo gas: a model in 
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which the hot gas is heated in situ by extraplanar SNe 
and is contain ed in isolated supernova remnants (SNRs; 
iSheltonI [20061 ) ■ and a model of an SN-driven ISM, one 
feature of w hich is the transfer of hot gas from the disk 
to the halo (jJoung fc Mac Lowll2006D . In addition, we 
use our observed halo parameters to estimate the X-ray 
luminosity of the halo, and compare it to the predic- 
tions of disk galaxy fo rmation mo dels (iToft et all 120021 : 
iRasmussen et al.ll2009l: iCrain et al.ll2O10ir 

Our XMM-Newton observations are a subset of those 
used in the survey of Henley & Shclton (2010, hereafter 
Paper I), who measured the SXRB O vii and O viii 
intensities from 590 archival XMM-Newton observations 
between I = 120° and I = 240°. The observations used 
here were chosen because they should be less affected 
by solar wind charge exchange (SWCX) emission 
(jCravensI l2000t ). w hich is a time-varyi ng conta minant 
of SXRB spectra (ICrav ens et al.l 120011: [Snowden~et alj 



20041: iFuiimoto et al.l~! 2007: Koutroumpa et a.1.1 



Kuntz fc SnowdenI I200& 



Henlev fc Sheltoa 120081) 



iCarter fc SembavT 
Although the 26 



2007 



2008 



XMM- 



Newton observations used here are only a small subset 
of the observations used in Paper I, this is a larger 
number of observations than has been used in pre- 
vious studies of the SXRB and the hot ISM with 
CCD-resolution spectra 

The remainder of this paper is organized as follows. In 
Section [2] we present the details of our observations and 
give an overview of the data reduction (see Paper I for 
more details). Section[3]contains our spectral analysis, in 
which we use CIE models to determine the spectrum of 
the halo emission. In Section H] we compare the results 
of our spectral analysis with the predictions of various 
physical models for the origin of the hot halo gas. In par- 
ticular, the disk galaxy formation model, the extraplanar 
SN model, and the SN-driven ISM model are presented 
in Sections 14.11 14.21 and 14.31 respectively. We discuss 
our results in Section [SJ and c onclude with a summary in 
Section[6] Throughout we use I Anders fc Grevessd (|19890 
abundances. 

2. OBSERVATIONS 

Our sample of observations is taken from Paper I, 
which presents O Vii and O VIII intensities extracted 
from 590 XMM-Newton observations between I = 120° 
and I — 240°. The observations used here were selected 
by applying various filters to minimize the contamina- 
tion from SWCX emission. In particular, to minimize 
geocoronal and near-Earth heliospheric SWCX contam- 
ination, we rejected the portions of our XMM-Newton 
data taken when the solar wind proton flu>0 exceeded 
2 X 10* cm~^ s^^, and to minimize heliospheric SWCX 
contamination we used only observations of high ecliptic 
latitudes (/3 > 20°) taken during solar minimum^ See 
Section 2 of Paper I for more details about SWCX and 
the filters we used to reduce its effects. 

After applying these filters, 43 observations remained 
(see Section 5.3.1 of Paper I). As we are interested in 

** The solar wind proton flux data were obtained from OMNIWeb 
(littp://omniweb.gsfc.nasa.gov ). 

'^ Following Paper I, we used observations made after 
00:OOUT on 2005 Jan 01. This data was estimated using 
sunspot dat a from the National Geophysical Data Center 
(pittp://www.ngdc.noaa.gov/stp/SOLAR/ftpsunspotnumber.html^. 
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Figure 1. All-sky Hammer-Aitoff projection in Galactic coordi- 
nates, centered on the Galactic Anticenter, showing the XMM- 
Newton pointing directions used in this study. The gray bands 
indicate regions of the sky that were excluded from our sample 
- the darker band shows |b| < 30°, and the lighter band shows 
\l3\ < 20°. The dashed line outlines the Eridanus Enhancement. 
The three observations toward this feature were also excluded. 

the Galactic halo here, we removed 14 more observations 
at low Galactic latitudes (|6| < 30°). The locations of 
the 29 remaining observations on the sky are shown in 
Figure [TJ Three of these observations are toward the 
Eri danus Enhancement, a large. X-ray-bright superbub- 
ble ([Burrows et al.lll993HS^owden et al.lll995D . We also 
removed these 3 observations. The details of the 26 ob- 
servations that remain in our sample are summarized in 
Table m 

The data reduction was carried out using SAS version 
T.O.CQ and the XMM-Newton Extend ed Source Analysis 
Software!!! (XMM-ESAS) vers ion 2 (jSnowden fc Kuntd 
120071: iKuntz fc Sn owdenl l2008D. We oiilv use d data from 
the EPIC-MOS cameras (jTurner et alll2001[) . as the ver- 
sion of XMM-ESAS that we used cannot calculate the 
parti cle background for EPIC-pn data (Str iider et al.l 
l2001h . The data reduction method is described in full 
in Paper I. Here we outline the main steps. 

We cleaned and filtered each dataset using the XMM- 
ESAS mos-filter script. This script runs the stan- 
dard emchain processing script, and then uses the XMM- 
ESAS clean-rel program to identify and remove times 
affected by soft-proton flaring. As noted above, we also 
removed times when the solar wind proton flux exceeded 
2 X 10* cm^^ s^^. The usable exposures for the two MOS 
cameras are shown in Columns 5 and 7 of Table !T| 

We detected and removed point sources with a 0.5- 
2.0 keV flux greater than 5 x 10"^''' erg cm ~^ s"""^ 



This i s the approximate flux level to which iChen et al.l 
(|1997f ) removed sources when they measured the 
spectrum of the extragalactic background to be 
10.5(^/1 keV)"i-4S photons cm'^ s'^ sr^^ keV"^ (we 



use this model in our spectral analysis; see Section 13. ip 
Some observations contained sources that were too bright 
or too extended to be removed by the automated source 
removal. We removed such sources by hand, by excluding 
a circular region centered on each source. 

We extracted SXRB spectra from the entire field of 
view, minus any sources that were removed, and minus 
any CCDs that were in window m ode or that exhib- 
ited the anomalous state described by IKuntz fc SnowdenI 

^ http://xmm2.esac.esa.int/sas/7.0.0/ 

''' http://heasarc.gsfc.nasa.gov/docs/xmm/xmmhp_xmmcsas.html 
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Table 1 

XMM-Newton Observation Details 













MOSl 




MOS2 




Obs. Id 


Start Date 


I 


b 


Exposure 


n 


Exposure 


n 








(deg) 


(deg) 


(ks) (arcmin^) 


(ks) 


(arcmin^) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


1 


0304070501 


2005-11-08 


124.223 


60.304 


12.0 


486 


12.0 


510 


2 


0305290201 


2005-07-02 


124.578 


-32.485 


15.1 


478 


16.9 


572 


3 


0401210601 


2006-10-10 


133.225 


42.419 


17.7 


310 


17.6 


394 


4 


0404220101 


2006-11-01 


135.974 


55.981 


13.0 


488 


14.1 


505 


5 


0400560301 


2006-11-17 


138.279 


68.853 


51.5 


380 


51.5 


403 


6 


0400570201 


2006-11-25 


142.370 


51.705 


23.0 


463 


22.9 


475 


7 


0406630201 


2007-04-12 


151.186 


48.245 


8.5 


339 


8.4 


431 


8 


0303260201 


2005-04-07 


151.607 


51.006 


44.4 


411 


44.0 


583 


9 


0306060201 


2005-11-13 


151.829 


70.103 


53.5 


410 


54.8 


511 


10 


0306060301 


2005-11-15 


151.831 


70.103 


15.3 


415 


15.6 


511 


11 


0303720601 


2005-04-25 


161.440 


54.439 


23.5 


382 


22.9 


398 


12 


0303720201 


2005-04-13 


161.441 


54.439 


25.8 


378 


26.2 


470 


13 


0200960101 


2005-03-28 


162.721 


41.656 


56.9 


453 


57.0 


465 


14 


0301340101 


2006-04-12 


167.648 


37.517 


12.8 


487 


13.0 


512 


15 


0306370601 


2005-04-24 


170.477 


53.178 


9.6 


496 


10.1 


583 


16 


0402780701 


2007-03-28 


171.132 


32.731 


14.3 


422 


14.6 


577 


17 


0304203401 


2006-06-11 


175.807 


63.353 


8.5 


371 


8.3 


391 


18 


0406610101 


2006-11-05 


179.356 


59.942 


10.4 


402 


10.8 


485 


19 


0400830301 


2006-10-30 


182.658 


42.566 


45.0 


493 


44.7 


511 


20 


0301651701 


2006-06-20 


197.309 


81.121 


12.3 


473 


12.2 


495 


21 


0300630301 


2006-01-19 


209.821 


-65.146 


14.4 


477 


14.5 


499 


22 


0312190601 


2006-01-28 


213.849 


-50.846 


11.3 


391 


11.2 


477 


23 


0301330401 


2006-02-13 


226.946 


-45.906 


19.5 


414 


19.6 


583 


24 


0312190701 


2006-01-28 


236.040 


-32.583 


11.1 


483 


10.8 


571 


25 


0302500101 


2005-08-09 


237.074 


-65.638 


21.9 


413 


23.5 


583 


26 


0307001401 


2006-02-13 


237.615 


-34.679 


7.8 


489 


8.3 


584 



Note. — The observations are in order of increasing I. Column 1 contains the XMM-Newton 
observation ID. Column 2 contains the observation start date, in YYYY-MM-DD format. 
Columns 3 and 4 contain the pointing direction in Galactic coordinates. Column 5 contains the 
usable MOSl exposure time that remains after the filtering mentioned in Section^ Column 6 
contains the solid angle, f2, from which the MOSl SXRB spectrum was extracted, after the 
removal of sources and unusable CCDs. Columns 7 and 8 contain the corresponding data for 
M0S2. 



(|2008f) . The solid angles from which each spectrum was 
extracted are shown in columns 6 and 8 of Table [1] We 
binned the spectra such that each bin contained at least 
25 counts. We created redistribution matrix files (RMFs) 
and ancillary response files (ARFs) using rmfgen and 
arf gen. 

We used the XMM-ESAS xmm-back program to calcu- 
late the quiescent particle background (QPB) spectrum 
for each observation. The QPB spectra were constructed 
from a database of filtcr-wheel-closed data, scaled using 
data from the une xposed corner pixels tha t lie outside 
the field of view (see lKuntz fc Snowdenlf2008l for more de- 
tails of the modeling of the QPB spectrum). The QPB 
spectra were subtracted from the corresponding SXRB 
spectra before we carried out our spectral analysis. 

3. SPECTRAL ANALYSIS 

3.1. Model Description 

We analyzed the spectra from each of our 26 XMM- 
Newton observations using a multicomponent model of 
the SXRB. The model consisted of the following compo- 
nents: a foreground emission component (representing 
LB emission and/or SWCX emission that remains in our 
spectra, despite our efforts to minimize this contamina- 
tion), a Galactic halo component, an extragalactic com- 
ponent, and instrumental background components. 

We used a .Raymond fc Smith (|1977. and updates) 



model with T '^ 10^ K to model the foreground emis- 
sion in our analysis, because such a model provides 
a good f it to data from the ROSAT All-Sky Sur- 
vey ( e.g., iSnowden etlllllQQSl 120001 : iKuntz fc SnowdenI 
[2000f) . We fixed the temperature (T = 1.2 x 10^ K) 
and normalizat i on of this component using data from the 
ISnowden etlUI (|2000( ) catalog of shadows in the ROSAT 
All-Sky Survey. Following Paper I, we calculated the 
foreground R12 (1/4 keV) count-rate for each XMM- 
Newton sightline by averaging the foreground count-rates 
of the 5 nearest shadows, weighted by their inverse dis- 
tances from the XMM-Newton sightline, i.e.. 



Foreground R12 count-rate : 






(1) 



where Ri is the foreground R12 count-rate for the ith 
shadow, and 9i is the angle between the center of the ith 
shadow and the XMM-Newton sightline. Over the whole 
set of XMM-Newton sightlines, the median value of 9i is 
6.2° (lower and upper quartiles: 4.1° and 7.7°, respec- 
tively). The median minimum and maximum values of 
6i for each sightline are 2.8° and 7.4°, respectively. Only 
three sightlines have Oi > 7° for ah i: obs. 0305290201 
[e, = 7.4°-9.3°), obs. 0301340101 {9, = 8.6°-14.5°), and 
obs. 0402780701 {6, = 11.2°-14.7°). We used the above- 
calculated count-rate to determine the normalization of 
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the foreground component for the sighthne in question 
(assuming T — 1.2 x 10^ K). This normahzation was 
held fixed during the subsequent spectral fitting. 

For the Galactic halo, which is the component of the 
SXRB that we are interested in here, we u sed a single- 
temperature (IT) Raymond &: Smith (|1977l and updates) 
model. While such a halo model is inadequate for ex- 

[ laining all the available X-ray and far-ultraviolet dat a 
Yao fc Wand[2007l: IShelton et al][2007l : ILei et al.l[2Q0l . 
it is adequate for characterizing the X-ray emission in 
the XMM-Newton band. The temperature and emission 
measure of this component w ere both free parame ters in 
the spectral fitting. We used a lRavmond fc Smithl model, 
instead of, say, an APEC model, because the codes that 
we used to calculate X-ray spectra fr om hydrodynamical 
models (see Section |4]) also use the iRavmond fc Smith! 
code. In our analysis, the temperature and emission mea- 
sure of the Galactic halo component were free to vary. 

We modeled the extragalactic background as a power- 
law with a photon index T = 1.46 and a normaliza- 
tion of 10.5 phot ons cm^^ s^^ sr^^ keV^ at 1 keV 
(jChen et al.l 119971 ) . The extragalactic and halo com- 
ponents were both subject to absorption. For each 
sightline, we used the HEAsoft nh tool to obtain the 
H I column density from the Leiden- Argentine-Bon n 
(LAB) Survey of Galactic H i (jKalberla et al.l [20051 ). 
We us ed photoelectric absorption cross-sections from 
iBaluci hska-Church fc McCa mmon (19 9 2), w ith an up- 
dated He cross-section from iYan et al.l (jl998f) . 

As well as the above SXRB components, we included 
components to model parts of the particle background. 
These components were independent for the two MOS 
detectors. The QPB spectrum includes two bright flu- 
orescent lines from aluminum and silicon at 1.49 and 
1.74 keV, respectively. The procedure for calculating the 
QPB spectrum mentioned in Section[2]cannot adequately 
remove these lines. Instead, the QPB spectrum was in- 
terpolated between 1.2 and 1.9 keV, and we added two 
Gaussians to our spectral model to model these lines. In 
addition, despite the data cleaning described in Section[2l 
some residual soft-proton contamination may remain in 
the data. We modeled this contamination using a power- 
law whi ch we did not fold through the instrumental re- 
sponse (|Snowden fc Kuntd [20071 boi O^. Following sug- 
gestions in the latest vers ion of the XMM-ESAS manual 
(jSnowden fc Kuiita '2010. which pertains to a later ver- 
sion of the software than the one we used), we placed con- 
straints on the spectral index of this power-law (specif- 
ically, soft limits at 0.5 and 1.0, and hard limits at 0.2 
and 1.3). 

Originally, as in Paper I, we used a broken power-law 
to model the soft protons . We fixed the break at 3.2 keV 
(jKuntz fc SnowdenI |2008[) . but we did not impose any 
special constraints on the spectral indices. We find that 
the halo temperatures are generally not significantly af- 
fected by our choice of soft proton model, but the emis- 
sion measures and surface brightnesses are typically 30% 
higher if we place constraints on the power-law spectral 
index. However, these differences do not affect the con- 
clusions of this paper. Throughout this paper, we use the 
newer set of spectral fit results, obtained using an unbro- 
ken power-law with constraints on the spectral index to 
model the soft protons. 



We carried out our spectral analysis using XSPEC0 
version 12.5.0. We analyzed each of our 26 observations 
individually, fitting the above-described model to the 
0.4-5.0 keV MOSl and M0S2 spectra simuhaneously. 

3.2. Systematic Errors 

Fixing the spectra of the foreground and extragalactic 
components may introduce systematic errors to our fit- 
ting (i.e., the true spectra of these components may differ 
from our assumed spectra, which may in turn affect the 
measured halo parameters). Here, we estimate the sizes 
of these systematic errors. 

In the case of the foreground model, we fixed the 
spectrum because otherwise there would be a degen- 
eracy between the foreground and background inten- 
sities. This deg eneracy can be o vercome b y shadow- 
ing observations (iSmith et al.ll2007l: [G aleazzi et al.ll2007l; 
Henlev fc SheltonI 120081: iGupta et ali i2009i: ILei erall 
20091), but such an analysis is not possible here. Here, 
we fixed the foreground normalization by extrapolating 
the foreground spectrum from the R12 band into the 
XMM-Newton band. However, as the relative contribu- 
tions of LB and SWCX emission are likely to differ in 
these bands, and these two emission mechanisms have 
different spectra, such an extrapolation may lead to an 
incorrect estimate of the foreground normalization. We 
estimated the size of this systematic effect by reanalyzing 
each sightline, using the median foreground R12 inten- 
sity (600 counts s~^ arcmin"^) to fix the normalization 
for every sightline. We used the median differences be- 
tween our original results and these new results as our 
estimate of the systematic errors due to the foreground 
normalization being fixed (we adopted the same system- 
atic errors for each sightline). The estimated systematic 
errors are ±0.08 x 10^ K for the halo temperature and 
±0.04 dex for the halo emission measure. 

We fixed the spectrum of the extragalactic background 
component because otherwise there would be a degen- 
eracy between this component and the power-law com- 
ponent used to model the soft protons. However, the 
normalization of this component may vary from field 
to field, due to statistical fiuctuations in the number 
of unresolved sources that comprise the extragalactic 
background. We estimated the size of the field-to- 
field variation in the extragalactic backgro und using the 
0.5-2. keV source flux distribution from iMoretti et al.l 
(|2003f ). Given that we removed sources with fluxes 
greater than 5 x 10^^* erg cm^^ s^"'^, we estimate that 
the 0.5-2.0 keV extragalactic surface brightness varies by 
roughly ±10% from field to field, assuming a typical field 
of view of 480 arcmin^. We therefore repeated our analy- 
sis for each sightline with assumed extragalactic normal- 
izations of 9.5 and 11.6 photons cm^^ s~^ sr^^ keV^ . 
We used the differences between the original results and 
these new results to estimate the systematic errors for 
each sightline due to our fixing the normalization of the 
extragalactic background. The systematic errors on the 
halo temperature and emission measure are typically less 
than ±0.2 x 10^ K and ±0.2 dex, respectively. 

An additional possible source of systematic error is 
the expected steepening of the extragalactic background 

* http://hcasarc.gsfc.nasa.gov/docs/xanadu/xspcc/ 
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below ^2 keV. In the 3-10 keV energy range, the ex - 
tragalactic photon index F « 1.4 (|Marshan et al.lll980l ). 
This photon index also pro vides good fits to S XRB spec- 
tra at lower energies (e.g.. f Chen et al.lll997L the source 
of our assumed extragalactic background model). How- 
ever, the summed spectrum of the individual sources that 
comprise the extr agalactic background is steeper below 
~2 keV (r = 1.96: lHasinger et al.|[l993l ). Failing to take 
this steepening into account would cause us to underesti- 
mate the low-energy extragalactic surface brightness and 
thus overestimate the halo surface brightness. We esti- 
mated the size of this systematic effect by using the dou- 
ble broken p ower- law model for the extragalactic back- 
ground from iSmith et al.l (|2007D . Both broken power- 
laws have F2 = 1.4 above the break energy of 1.2 keV. 
The first component has Fi = 1.54 below 1.2 keV, and 
a normalization of 5.70 photons cm~^ s""*^ sr"-*^ keV~ . 
The second component has Fi = 1.96 and a nominal nor- 
malization of 4.90 photons cm~^ s^^ sr^^ keV^^. Note 
that the sec ond normalization was a free parameter in 
ISmith et al.l 's (2007) analysis, and they obtained a value 
roughly 50% lower. However, because of the aforemen- 
tioned degeneracy between the extragalactic background 
and the soft proton component, we fix the normaliza- 
tion at is nominal value. Note also that the spectrum of 
the extragalactic background depends on the flux level to 
which sources are rem oved; it is not clear to what source 
removal threshold the ISmith et al.l ((2007 ) model is ap- 
plicable. Nevertheless, we proceeded by repeating our 
analysis for each sightline, using this new extragalactic 
background model. As with the normalization of the ex- 
tragalactic background, we used the differences between 
the original results and these new results to estimate 
the systematic errors for each sightline. The systematic 
errors on the halo temperature are typically less than 
±0.1 X 10^ K. As expected, this new extragalactic model 
generally yields lower halo emission measures, although 
the difference is typically less than 0.2 dex. 

To calculate the systematic errors on the halo param- 
eters for each sightline, we added the three systematic 
errors discussed above in quadrature. We report these 
combined systematic errors alongside the statistical er- 
rors in the following section. 

3.3. Spectral Fit Results 

The results of our spectral analysis are presented in Ta- 
ble [2] In particular, columns 4 and 5 give the best- fitting 
halo temperature and emission measure (E.M.) for each 
observation, and column 7 gives the intrinsic 0.4-2.0 keV 
surface brightness of the halo component. For the tem- 
perature and the emission measure, we present both the 
Icr statistical errors and the estimated systematic errors 
discussed in the previous section. 

The best-fit halo parameters are plotted against Galac- 
tic latitude in Figure [5] The error bars show the statisti- 
cal and systematic errors added in quadrature. The ob- 
servations with large temperature error bars at |&| w 33° 
and 35° are 0402780701 and 0307001401 (numbers 16 
and 26, respectively, in Tables [T] and [2|) . Note that the 
halo is faint in these directions (indeed, these observa- 
tions yield the two lowest 0.4-2.0 keV surface bright- 
nesses), and that the exposure times are not unusu- 
ally long. It is therefore unsurprising that the halo 



temperatures are poorly constrained for these observa- 
tion. The other temperatures are generally well con- 
strained; they are typically ^^(1.8-2.4) x 10^ K, and 
are fairly constant across the sky. The halo emission 
measures mostly lie in the range ^0.0005-0.006 cm^^ 
pc, although there is significant variation within that 
range. Correspondingly, there is also significant varia- 
tion in the intrinsic X-ray surface brightness, within the 
range ^(0.5-5) x 10^^^ erg cm~^ s~^ deg^^. 

One observation, 0406630201 (number 7 in Table [5]) , 
gives a significantly higher halo temperature than the 
other observations (11.2 x 10^ versus ~2 x 10^ K). 
The M0S2 spectrum and best-fit model for this ob- 
servation are shown in Figure [31 For comparison, Fig- 
ure [3] also shows the spectrum from a nearby observation 
(0303260201; number 8 in Table [5]) that yields a more 
typical halo temperature (1.8 x 10^ K). The halo com- 
ponent for obs. 0406630201 is faint, and the temperature 
may be less well constrained than the formally calculated 
error bar implies. We tried repeating our analysis of this 
observation with the halo temperature fixed at 2 x 10^ K 
(i.e., similar to the temperatures found from most of the 
other observations). The best-fit emission measure for 
this new halo component was essentially zero (3ct up- 
per limit: 0.0022 cm~^ pc). This new model yielded 
an acceptable fit: x^ ~ 309.59 for 287 degrees of free- 
dom. Therefore, although the formal best-fitting tem- 
perature is 11.2 X 10^ K, this spectrum is also consistent 
with a ^2 X 10^ K halo with a small emission measure 
(<0.002 cm~^ pc). The upper limit on the emission mea- 
sure is not unusually small compared with some of the 
other observations. 

Although there is this one anomalous temperature 
among our 26 XMM-Newton observations, it should not 
affect our subsequent analysis. In Section JU we will com- 
pare our observations of the halo with the predictions of 
various physical models. Rather than looking at individ- 
ual sightlines, we look at the whole population of observa- 
tional results, and compare histograms of observational 
properties with corresponding histograms derived from 
the models. Therefore, as long as the majority of our 
observations yield reasonably accurate halo properties, a 
single outlying anomalous result should not significantly 
affect the comparison of the observations with the mod- 
els. 

4. COMPARING THE OBSERVED HALO X-RAY SPECTRA 
WITH HYDRODYNAMICAL MODELS 

In this section, we compare the X-ray spectral prop- 
erties of the halo inferred from our XMM-Newton ob- 
servations with those predicted by various hydrodynam- 
ical models. In particular, in Section 14.1! we examine a 
disk galaxy formation model, in which extragalactic ma- 
terial is heated as it falls into the Galaxy 's potential well 
(jToft et al.ll2002t IRasmussen et al.ll2009h . In Section ill 
we examine a model in which the hot halo gas is heated in 
situ by isolated extraplanar SNe ()Sheltonll2006[) . In Sec- 
tion |T31 we examine a model in which th e ISM is heated 
and stirred by multiple SNe (Joung & Mac Low! [2006'). 
Unlike the previous model, the SNRs are not assumed to 
evolve in isolation, and the model includes the movement 
of hot gas from the disk into the halo. 

4.1. Disk Galaxy Formation Model 
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Table 2 

Spectral Fit Results 



Obs. ID Foreground R12 rate A^h Halo T 

{ROSAT units) (lO^o cm-2) (lO" K) 
(1) (2) (3) (4) 



Halo E.M. 
(10-3 cm-6 pc) 

(5) 



X /dof 


50.4-2.0 




(10-12 erg cm-2 s 


(6) 


(7) 


405.97/362 


1.66 


475.90/440 


3.73 


452.03/395 


1.75 


469.00/422 


1.40 


693.09/594 


2.71 


504.83/529 


1.57 


285.43/286 


1.89 


607.54/593 


3.05 


587.85/594 


2.13 


405.23/412 


2.60 


498.38/445 


1.93 


494.35/494 


2.68 


606.12/594 


2.01 


349.15/403 


2.73 


353.40/340 


2.11 


400.90/451 


0.58 


271.14/239 


3.91 


367.42/368 


3.25 


630.86/594 


1.80 


386.43/376 


4.65 


380.59/413 


0.84 


322.29/326 


1.84 


457.39/496 


1.54 


375.97/363 


3.21 


446.89/476 


5.48 


380.72/335 


0.51 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 



0304070501 
0305290201 
0401210601 
0404220101 
0400560301 
0400570201 
0406630201 
0303260201 
0306060201 
0306060301 
0303720601 
0303720201 
0200960101 
0301340101 
0306370601 
0402780701 
0304203401 
0406610101 
0400830301 
0301651701 
0300630301 
0312190601 
0301330401 
0312190701 
0302500101 
0307001401 



743 
293 
593 
606 
643 
586 
577 
555 
687 
687 
641 
641 
615 
575 
742 
471 
763 
806 
527 
542 
698 
427 
417 
573 
597 
666 



0.790 
5.83 
3.45 
1.72 
1.60 
0.642 
1.00 
0.640 
1.26 
1.26 
1.21 
1.21 
1.94 
3.31 
0.656 
4.54 
1.19 
1.65 
1.74 
1.73 
2.13 
2.29 
2.73 
1.75 
3.01 
2.63 



r, Q-,+0.08 +0.16 
^••JJ^_0.24 -0.24 
r, 1Q+0.07 +0.21 
^■^^='-0.08 -0.16 
r, 9C- + 0.11 +0.16 
^■^3-0.23 -0.21 
, QC-+0.20 +0.72 
l-»O-0.19 -0.17 

2.10t°°6 ±0.16 

r, r,O + d.09 +0.17 

^•'JO-O.ia -0.24 
11 1S+1-04 +0.16 

11.10_Q gQ _]^ 41 

1 OO+0.02 +0.16 
1-OJ-0.06 -0.17 
1 /ffi+0-03 +0-49 
l-^D-0.04 -0.16 
9 nK+0-08 +0.16 
^•I-ID-O.IS -0.21 
r, -,1+0.12 +0.32 
^•11-0.09 -0.17 
r, ns+O-OS +0.17 
^■U0_o.o6 -0.16 

2.20tnSa ±0.17 



1.63l° 



1.99 



0.08 
,+0.06 



±0.16 



-0.11 
-I 79+0.12 +0.16 
^•'^-0.10 -0.18 

r, r|O + 0.37 +9.89 

^•UO_o 4g _o.20 

1 Qfi+0-09 +0.17 
l-yO-O.lO -0.16 

14+0.42 +0.61 
■J-l''-0.23 -0.22 

1 Q1+0.10 +0.22 
J^-^J^-0.06 -0.18 

1.83 ±0.05 ±0.16 

1 QQ+0.26+0.17 
l-»y-0.36 -0.19 

2.28t°" ±0.19 

r, 1Q + d.ll +0.41 

^•iy-0.16 -0.16 
9 nA+0.07 +0.16 
^•'-'^-0.10 -0.19 
r, qc+O.OS +0.19 
^•■J3-0.06 -0.16 
n O1+0.90 +7.62 
■^•■^-'--1.09 -0.47 



3.96 



i.solg 



1.78 



1.79Z 
1.641 



5.541 



3.051° 

1+0 

4.49lf 

3.0llg 
2.161° 
3.07lg 
2.121° 
3.35l° 
3.55l° 



0.69^ 



4.97 



+0 



-0 

2.39lg 
6.79lg 

+0 
-0 

5.27l° 



2.45^ 



1.031 
1.851 
1.63 
3.78^ 



0.37^ 



.20 +0.29 
25 -0.64 
27 +0.78 
.45 -0.76 
.27 +0.45 
34 -0.52 

18 +0.32 
.52 -1.35 
.17 +0.56 

19 -0.64 
.23 +0.40 
.19 -0.53 
.16 +0.29 

49 -0.56 

33 ± 1-01 
.27 +0.98 
52 -3.35 

16 +0.54 
.49 -0.57 
.26 +1.09 

25 -0.97 

22 +0.88 
.36 -0.53 
.11 +0.39 

19 -0.45 
.29 +0.74 
.42 -0.64 
.35 +0.81 
0.58 
.37 +0.56 
.36 -0.45 
.47 +0.88 

54 -1.16 
.23 +0.44 
.69 -0.92 
.28 +0.79 

17 -0.78 
.37 +1.20 
.46 -1.39 
.43 +0.58 

32 -0.32 
.25 +0.73 
.39 -0.72 
.25 +0.36 

24 -0.84 
.31 +1.04 
.33 -0.85 
.13 +0.94 

36 -1.50 
.19 +0.07 
.23 -0.29 



Note. — Column 1 contains the XMM-Newton observation ID. Column 2 contains the foreground R12 (1/4 kcV) count-rate in ROSAT units 
(IQ- counts s- arcmin— ). This count-rate was derived from the lSnowden et all ( 120001) catalog of SXRB shadows, and was used to fix the normal- 
ization of the 1.2 X 10 -K foreground component. Column 3 contains the H I column density jKalberla et aUI2005l) that was used to attenuate the halo 
and extragalactic components. Columns 4 and 5 contain the best-fit halo parameters (E.M. — Jn^dl). In each case, the first error indicates the la 
statistical error, and the second error indicates the estimated systematic error due to our assumed foreground and extragalactic spectra (sec Section l3.2l 
for details). Column 6 contains x ^^^^ th® number of degrees of freedom. Column 7 contains the intrinsic 0.4—2.0 keV surface brightness of the halo 
component. 



Cosmological smoothed particle hydrodynamics (SPH) 
simulations of disk galaxy formation predict that such 
galaxies should be surrounded by ext ended hot halos 
(r ~ 10s of kpc T ^ few x 10^ K; iToft et al.l [20021: 
iRasmussen et al.]|2009[ ) . These halos con tain a significant 
fracti on of the galactic baryonic mass (jSommer-Larsenl 
|2006[) . In this subsection we compare the predictions of 
such models with our XMM-Newton observations. 

J. Rasmussen (2009, private communication) has 
kindly provided us with 0.3-2.0 keV luminosities and 
emission-weighted mean temperatures derived from 
the SPH galaxy formation simulations described in 
IRasmussen et al.l (|2009( ). IRasmussen et al.l ([2009^ point 
out that the X-ray emission from the hot gas particles 
can be artificially boosted by nearby cold, dense gas 
particles (within an SPH smoothing length, /isph)- As 
/igPH w 1.5 kpc in these simulations, the emission from 
the disk can be particularly adversely affected. There- 
fore, a cylindrical region around each galactic disk within 
|2;| = 2 kpc and r = 15 kpc was excluded from the cal- 
culation of the X-ray properties. For each galaxy, two 
sets of values were extracted: one extracted from within 



a spherical region of radius 40 kpc (t hese are the v alues 
shown in Figures 5 and 6 in .Rasmussen et al.ll2009l ). and 
one extracted from within a (100 kpc)^ box. The X-ray 
luminosities and temperatures for each model galaxy are 
shown by the triangles in Figure |4l 

In order to compare the model predictions with our 
observations, we must first derive a luminosity for the 
Galactic halo from our observations. To do this, we 
must assume some geometry for the halo emission. In 
the following, we consider both spherical and cylindrical 
geometries. Of the two, a spherical geometry is probably 
the more appropriate for comparison with the extended 
hot halo predicted by the model, especially as a cylindri- 
cal region around each model galactic disk was excluded 
before the X-ray properties were extracted. 

For the spherical halo geometry, we assume that the 
halo emission comes from a uniform sphere of radius 
i?sph- If the intrinsic surface brightness, ^x, along some 
sightline is typical for the whole galaxy, the luminosity 
per unit volume of the halo is 47rS'x/A, where A is the 
path length through the spherical halo. This path length 
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Figure 2. Halo (a,b) temperature, (c) emission measure, and (d) intrinsic 0.4-2.0 keV surface brightness plotted against Galactic latitude. 
Panel (b) shows the same data as panel (a), but with a narrower y-axis range. The error bars show the statistical and systematic errors 
added in quadrature. The errors on the surface brightness are derived from the errors on the emission measure. 
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Figure 3. MOS2 spectra from (a) the observation with an exceptionally high halo temperature of 11.2 X 10 K, obs. 0406630201, and (b) 
a more normal observation, obs. 0303260201, showing the best-fitting model (in red) and the individual model components. Components 
of the particle background are plotted with dashed lines. 



is a function of the viewing direction: 



X = Rq cos lcosb+ JR'^pi^ - (1 - cos2 1 cos^ 6)i?|, (2) 

where Rq = 8.5 kpc is the radius of the solar circle. The 
X-ray luminosity is then given by 



Lx = 



IGtt^ Rsph^^ 



(3) 



3 A 

For the cylindrical halo geometry, if the intrinsic sur 



face brightness along some sightline is typical for the 
whole galaxy, the luminosity per unit area of the halo 
integrated over the vertical direction is 2 x 47rsin(|6|)S'x- 
The initial factor of 2 takes into account the halo above 
and below the disk. If we then assume that the halo 
emission originates within a cylindrical region of radius 
i?cyi, the total X-ray luminosity of the Galactic halo is 

Lx^8n^smi\b\)Rl^iSx. (4) 

We have calculated intrinsic 0.3-2.0 keV surface bright- 
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Figure 4. Halo temperatures and 0.3-2.0 keV luminosities de- 
rived from our XMM-Ne wton observations and from disk galaxy 
formation simulations iRasmussen et al.l2009 ) . The triangles show 
the emission-weighted mean temperatures and halo luminosities 
predicted by the model. These datapoints are labeled with the 
model galaxies' circular velocities, Vc, in km s~^. For each model 
galaxy, the lower datapoint shows the values extracted from within 
a spherical region of radius 40 kpc, and the upper datapoint shows 
the values extracted from within a (100 kpc)'^ box. In both cases, a 
cylindrical region around the galactic disk was excluded. The solid 
circle and the diagonal cross show the median observed values. The 
solid circle denotes the median halo luminosity derived from the ob- 
servations assuming a spherical emission geometry (Equation lO); 
the errorbars indicate the lower and upper quartiles. The diago- 
nal cross denotes the median halo luminosity derived assuming a 
cylindrical emission geometry (Equation (O). 



nesses for the hafo from our best-fitting spectral models 
in Table El We then converted each surface brightness to 
a 0.3-2.0 keV luminosity, using Equations ^ and ^. In 
both cases we assumed emission radii (i?sph or i?cyi) of 
15 kpc. There is a large amount of scatter in the derived 
luminosities, spanning about an order of magnitude. The 
medians of these values are 1.9 x 10^^ erg s~^ (spherical 
geometry) and 1.2 x 10^^ erg s~^ (cylindrical geometry); 
these are our best estimates of the 0.3-2.0 keV luminos- 
ity of the Milky Way halo. However, it should be noted 
that our sightlines are all in directions away from the 
Galactic Center, and so if the halo emissivity increases 
toward the Galactic Center, the above values will under- 
estimate the Galactic luminosity. In addition, the lumi- 
nosity inferred assuming a spherical geometry will be an 
underestimate if the halo is more extended than our as- 
sumed i?sph — 15 kpc (note from Figure|3]that the model 
predicts that a significant fraction of the halo emission 
originates from r > 40 kpc). For comparison, analysis 
of ROSAT All-Sky Survey data has yielded halo X-ray 
luminosities of 7 x 10^^ erg s~^ in the 0.1-2.0 keV band 



(assuming that the quoted valu e applies to the whole 
ROSAT band; IPietz et all [TOOl and 3 x lO^^ erg s'^ 
in the 0.5-2.0 keV band ()Wanelll998[) . Assuming a halo 
temperature of 2 x 10^ K, these values correspond to 
0.3-2.0 keV luminosities of 4 x 10^^ and 5 x 10^^ erg s^^, 
respectively. 

Figure |4] compares the median observed halo tempera- 
ture and the above-mentioned median halo luminosities 
with the predictions of the disk galaxy formation simula- 
tions. For model galaxies similar in size to the Milky 
Way {vc « 220 km s~^), the predicted temperatures 
are in good agreement with the observations. However, 
the IRasmussen et al.l ()2009D model underpredicts the ob- 
served halo luminosity by at least an order of magnitude. 
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Figure 5. Density and temperature profiles for an SNR model 
with _Bo = 0.5 X 10^^ erg, Befi = 2.5 /xG evolving in an ambient 
medium with no = 0.01 cm^^. Note that the cool shell has formed 
by 250,000 yr. 

We will discuss these results in Section r5.2.1l 

4.2. Extra-planar Supernova Explosions 

Here, we consider a model in which the hot halo gas 
is heated locally by SNe above the Galactic disk. In this 
scenario, the observed hot gas is within isolated SNRs at 
a variety of heights above the disk and of a variety of ages. 
Superbubbles blown by clustered SNe and hot gas that 
has risen from the disk a,re excluded in this model. This 
model was developed bv iSheltoiil ()2006D . who found that 
it could explain a significant fraction of the high-latitude 
1/4 keV halo emission (excluding anomalously bright fea- 
tures such as the North Polar Spur). Here, we compare 
the model to the higher-energy emission observed with 
XMM-Newton. 

4.2.1. Model Description 

The model spectra discussed here were generated from 
one-dimensional Lagrangian hydrodynamical simulations 
of SNRs evolving at a variety of heights above the Galac- 
tic disk, and hence in a variety of a mbient d e nsitie s . The 
sim ulations are described fully in iSheltonI ()1998l 119991 
120061) ■ and include radiative cooling, thermal conduction, 
and an effective ambient magnetic field, Bcff, which ex- 
erts a non-thermal pressure, in addition to the ambient 
gas pressure. The ionization evolution in the shocked gas 
is modeled self-consistently. 

We will first concentrate on the models from ISheltonI 
(2006) with SN explosion energy Eq — 0.5 x 10^^ erg and 
-Boff = 2.5 /iG (corresponding to a non-thermal pressure 
Pnt = 1800 cm~'^ K), and then discuss varying these pa- 
rameters. We consider the models evolving in ambient 
densities no = 0.2, 0.1, 0.05, 0.02, 0.01, and 0.005 cm'^, 
corresponding to heights z = 190, 310, 480, 850, 1300, 
and 180 pc, u sing the interstellar density model from 

iFerrierd (11998 ). Note that we are ignoring the model 

from ISheltonI ( 2006[ ) at 2 = 76 pc, as we do not consider 
this to be in the halo. Density and temperature profiles 
from the model with no = 0.01 cm~'^ are shown in Fig- 
ure [5l The SN explosion blows a hot rarefied bubble in 
the ambient medium; this hot bubble is the source of 
the X-rays. T he mode l show n in Figure [S] corresponds 
to model A in ISheltonI ()1999D ; see that paper for more 
details. 

For each epoch of each SNR model, we calculated spec- 
tra for a range of impact parameters through the model 
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remnant, using the iRavmond &: Smitbl (|1977[ ) spectral 
code (updated by J. C. Raymond & B. W. Smith, 1993, 
private communication with R. J. Edgar). The spac- 
ing between the impact parameters, A6, was adjusted 
according to the size of the remnant. The spectral cal- 
culation takes into account the possible non-equilibrium 
ionization calculated during the hydrodynamical simula- 
tion. 

We wish to compare the distributions of observed tem- 
peratures and emission measures with those predicted 
by the model. In order to do this, we need to calculate 
the probability of a sightline passing through a remnant 
at a given height, of a given age, and at a given im- 
pact parameter. For example, a sightline is more likely 
to pass through a remnant at a larger impact parame- 
ter than a small impact parameter, and more likely to 
pass through a remnant closer to the disk, where the SN 
rate is larger. If R{zi,Z2) is the rate per unit area of 
SNe at heights between zi and Z2, then the probability 
P{zi, Z2,t, At,b, Ab) of intercepting a remnant between 
these heights with an age between t and t + At at an 
impact parameter between b and b + Ab is 

P{zi,Z2,t, At, b, Ab) = 2TTbR{zi,Z2)AbAt. (5) 

In order to cal culate the above probabilities, we followed 
iSheltonI (J2006I ) and divided the halo into 6 plane-parallel 
slabs. We assumed that all the SNRs within a given 
slab experience a single ambient density, and so can be 
represented by one of our 6 SNR models. We calculated 
i?(zi, Z2) for each slab by integrat ing the volume tric SN 
rate at the solar circle, r{z), from iFerrierd ()1998D : 

r{z)^n^{z) + ni[z) 

= (4.0e-l"l/325 P= + 14e-l"l/266 p<=) kpc'^ Myr-\ (6) 

where ri^{z) and rii{z) are the rates of type la SNe and 
isolated core-collapse SNe, respectively. The slab bound- 
aries were placed at the midpoints between the nomi- 
nal heights of the model remnants, except for the lower 
boundary of the lowest slab, which was placed at 130 pc 
(the scale height of the Galactic H i layer) , and the upper 
boundary of the highest slab, which was placed at infin- 
ity. The boundaries, ambient densities, and SN rates for 
the 6 slabs are shown in Table [S We calculated At from 
the timestamps of the output files from the hydrody- 
namical simulations. As noted above, Ab is the spacing 
between the impact parameters for which we calculated 
spectra from a given remnant. 

We used a Monte Carlo method to construct model 
spectra for 2000 sightlines, using the above-calculated 
probabilities. In our Monte Carlo simulation, 64% of the 
model sightlines intercepted no remnants, and 9% of the 
sightlines intercepted more than 1 remnant. For model 
sightlines intercepting more than 1 remnant, we summed 
the spectra of the individual remnants. 

4.2.2. Characterizing the SNR Spectra with IT Models 

In our analysis of the XMM-Newton observations, 
we modeled the halo X-ray emission with a single- 
temperature CIE plasma model, whereas the true halo 
emission is likely from plasma at a range of temperatures 
and in a range of ionization states. Similarly, the emis- 
sion predicted by this extraplanar SNR model is from 
plasma at a range of temperatures and in a range of 



ionization states. Therefore, in order to compare the 
predictions of this model with our observational results, 
we first characterize the model SNR spectra calculated 
above with IT CIE plasma models, by simulating XMM- 
Newton observations of the SXRB. 

Our procedure for characterizing the spectrum for 
each model sightline is as follows. We first multiplied 
the spectrum by a renormalization factor, km, in or- 
der to give a 0.4-2.0 keV surface brightness of 2.06 x 
10~^^ erg cm~^ s~^ deg~ . This value is the median 
intrinsic halo surface brightness inferred from the best- 
fit models in Table [2j We then subjected the model to 
an absorbing column TVh = 1.7 x 10^° cm~^; again, this 
is the median value used in our XMM-Newton analy- 
sis. To the absorbed SNR spectrum, we added compo- 
nents representing the foreground emission, the extra- 
galactic background, the soft protons, and the instru- 
mental lines. The parameters for the foreground emis- 
sion and the extragalactic background were taken from 
our observational analysis (see Section 13. ip . The nor- 
malization of the foreground component was chosen to 
give an R12 count-rate of 600 counts s~^ arcmin"^ (the 
median value from Table [2]). The parameters for the soft 
proton model and the instrumental lines were taken from 
the best-fit model for obs. 0301330401; this observation 
was chosen because it has close to the median level of 
soft proton contamination, as judged by the ratio of the 
observed 2-5 keV fiux to that expected from the extra- 
galactic background (F^^itai/ ^cxglv '^ee Paper I). 

We simulated an observation of the model spectrum 
by folding the above-described multicomponent model 
through the XMM-Newton response, assuming a typi- 
cal field of view of 480 arcmin^, and adding Poissonian 
noise to the spectrum, assuming a typical observing time 
of 15 ks. Our simulations also took into account the 
QPB spectrum. For each model sightline, we simulated 
a MOSl and a M0S2 spectrum. We binned the resulting 
spectra such that there were at least 25 counts per bin. 
Note that our assumed field of view is not as large as the 
full XMM-Newton field of view, as we removed bright 
sources from the fields of our observations, and for some 
observations not all chips were usable. 

We fitted the resulting simulated spectra with the same 
multicomponent model that we used in our observational 
analysis. In partic ular, the halo emis sion was modeled 
with an absorbed IRavmond fc SmithI model, with JVh 
fixed at 1.7 X 10^" cm~^. For each model sightline we 
fitted to the simulated MOSl and M0S2 spectra simulta- 
neously. Figure [6] shows a simulated XMM-Newton spec- 
trum from one of our model sightlines, along with the 
best-fitting model. Also shown is the IT CIE halo com- 
ponent of the best-fitting model, as well as the input SNR 
model. The input SNR spectrum is very different from 
that of a CIE plasma. The SNR plasma is relatively cool 
(as shown by the differential emission measure, plotted as 
the solid line in Figure [7|) and overionized, and the spec- 
trum exhibits strong recombination edges. If the SNR 
plasma were in CIE, instead of recombining, its emission 
would be 4 orders of magnitude fainter in the XMM- 
Newton band. However, when we use a IT CIE plasma 
to model the halo component in our simulated spectra, 
the fits are generally good. This is because much of the 
spectral detail in the input SNR spectrum is lost when 



10 



HENLEY ET AL. 



Table 3 

Densities and Supernova Rates in the Halo as a Function of Height 

Slab Height range Ambient density, no Nominal SNR height SN rate, R 

(pc) (cm~'^) (pc) (kpc~^ Myr"-') 



1 


130-250 


2 


250-395 


3 


395-665 


4 


665-1075 


5 


1075-1550 


6 


1550-00 



0.2 
0.1 
0.05 
0.02 
0.01 
0.005 



190 


1.1 


310 


0.83 


480 


0.76 


850 


0.36 


1300 


0.091 


1800 


0.022 




Energy (keV) 

Figure 6. Simulated XMM-Newton MOSl spectrum of the 
SXRB, created using the extraplanar SNR model for the halo emis- 
sion. The top panel shows the input SNR spectrum. The middle 
panel shows the simulated SXRB spectrum (crosses). The upper 
solid line shows the best-fit SXRB model, while the lower solid 
line shows the IT halo component from that model (the other 
model components have been omitted for clarity; cf. Figure [Sjl. 
The dotted line shows the input SNR spectrum folded through the 
XMM-Newton response. The bottom panel shows the residuals. 
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Figure 7. Differential emission measures for the input SNR 
model plotted in Figure [6] (solid line) and the input SN-driven 
ISM model plotted in Figure [TOl (dashed line). 



it is combined with the other components of the SXRB 
and folded through the XMM-Newton response. There 
is some excess emission at ^0.75 keV in the simulated 



spectrum in Figure |51 possibly due to O 



+ 7 



-^ 



0+6 



combinations (the recombination edge is at 0.74 keV), 
but it would not be easy to unambiguously identify such 
a feature in an observed spectrum as recombination emis- 
sion. 

We used the resulting fit parameters to compare with 
the observations. The temperatures were taken directly 
from the fits, while the best- fit emission measures were 
first divided by the relevant values of km before compar- 
ing with the observations. In what follows, we refer to 
these values as "X-ray temperatures" and "X-ray emis- 
sion measures" , to emphasize that they are derived from 
simulated X-ray observations, rather than being derived 
directly from the hydrodynamical data. 

4.2.3. Comparing the Extraplanar SNR Model with 
Observations 

Our model X-ray emission measures were derived for 
sightlines looking straight upward from the disk, i.e., to- 
ward |6| = 90°. Our observations, however, are toward 
a range of Galactic latitudes, and so will sample differ- 
ent amounts of halo material. We therefore multiply our 
observed emission measures by sin|6| before comparing 
them with the model predictions. This transformation 
assumes that the halo is, in a statistical sense at least, 
uniform in directions parallel to the disk. 

FiguresIHK and[8b show histograms comparing the halo 
X-ray temperatures and emission measures predicted 
by the extraplanar SNR model with the corresponding 
halo properties obtained from our XMM-Newton ob- 
servations. (Note that obs. 0406630201, which yielded 
T = 11.2 X 10^ K, is not shown in FigurelDi.) 

The X-ray temperatures predicted by the model are in 
reasonably good agreement with the observed tempera- 
tures. However, the emission predicted by this model is 
clearly too faint: the predicted X-ray emission measures 
are typically an order of magnitude smaller than the ob- 
served values. Among the model sightlines that intercept 
extraplanar SNRs, the median X-ray emission measure is 



1.37 



+0.09 



dex smaller than the median observed value of 



-0.46) X 10" 



-0.12 

E.M. sin |6|: 0.081 x lO^^ versus (1.90^°'^^^ 
pclj As the predicted X-ray temperatures are in rea- 
sonable agreement with the observed temperatures, the 
emission measure result implies that the X-ray surface 



^ Here and in Section [4.3.2l the errors indicate the 90% bootstrap 
confidence interval on the observed median. 
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(a) Eo^O.SxIo'^' erg, Betf = 2.5|iG 
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(b) E„ = 0.5 X 1 0" erg, B,„ = 2.5 jiG 



Model I 



Observed 
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'1 
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Figure 8. Histograms comparing the halo X-ray temperatures 
(panel a) and X-ray emission measures (panels b and c) predicted 
by the extraplanar SNR model (solid gray) with those obtained 
from our XMM-Newton observations (black outline). The observed 
emission measures have been multiplied by sin |6| (see text for de- 
tails). The model temperatures and emission measures were de- 
rived b y cha racterizing the model spectra with IT models (see 
Section l4.2.2l l. Panels (a) and (b) show model values derived from 
simulations with Eq = 0.5 X 10^^ erg and B^g = 2.5 /iG, with 
a Galactic SN rate given by Equation ^. Panel (c) shows model 
emission measures derived from simulations with Eg = 1 X 10^^ erg 
and Bcff = 5.0 ^G, with a Galactic SN rate twice as large as that 
given by Equation ^. The model histograms have been rescaled 
so they cover the same area as the corresponding observational 
histograms. Note that 64% of the model sightlines in panels (a) 
and (b), and 42% of the model sightlines in panel (c) intercept no 
model remnants. Such sightlines have an undefined X-ray temper- 
ature and zero X-ray emission measure, and thus do not appear in 
the above plots. 



Table 4 

X-ray Emission Measures Predicted by Different 
Extraplanar SNR Models 



Eo B,g 

(10^1 erg) (ajG) 



SN rate'' 



Median E.M. 



(10- 



pc) 



0.5 





Ix 


0.028 


0.5 


2.5 


Ix 


0.081'' 


0.5 


5.0 


Ix 


0.19 


1.0 


5.0 


Ix 


0.21 


1.0 


5.0 


2x 


0.23<= 


Observed 






1 qn+0-4ld 

t-yu_o.46 



Note. — The emission measures were obtained by 
fitting to the spectra above 0.4 keV. 
^ Relative to the rate given by Equation JGl . 
^ Model used in Figures [SJi and|8j5. 
'^ Model used in FigurejSJ;. 

'^ Median of E.M.sin|b|, with 90% bootstrap confi- 
denee interval. 

brightnesses are also typically underpredicted by an or- 
der of magnitude. Note also that our Monte Carlo sim- 
ulation predicts that about two-thirds of the sightlines 
would intercept no extraplanar SNRs, whereas we ob- 
serve hot halo gas on most, if not all, of our sightlines. 

The model values in Figures [8^ and |8}d were calcu- 
lated from the SNR simulations in Shelton (2006) with 
Eo ^ 0.5 X 10^1 erg and B^e = 2.5 ^iG. We find that 
increasing Eq, i?eff, or the assumed SN rate all increase 
the predicted X-ray emission measures. Table |4] shows 
the median X-ray emission measures predicted by ex- 
traplanar SNR models with different explosion energies, 
ambient magnetic fields, and SN rates. As can be seen, 
increasing BcS has the largest effect on the predicted X- 
ray emission measures. However, it should be noted that 
the median model values in Table |4] are only for the sub- 
set of sightlines that intercept at least one SNR. Increas- 
ing Eq or the SN rate increases the fraction of sightlines 
that intercept at least one SNR, while increasing Bcs 
decreases that fraction. 

Figure |8}: shows the histogram of X- ray emission mea- 
sures predicted by the iSheltonI (|2006[ ) simulations with 
Eo = lx 10^1 erg and B^s = 5.0 //G, with a Galactic SN 
rate that is twice that given by Equation ([6]) (i.e., the 
5th model in Table |4]). The histogram of observed values 
is also plotted for comparison. The model still signifi- 
cantly underpredicts the observed emission measures (it 
underpredicts the median value by 0.911q'°2 dex). Fur- 
thermore, our new Monte Carlo simulation predicts that 
^40% of the sightlines would intercept no extraplanar 
SNRs - we reiterate that we observe hot halo gas on 
most, if not all, of our sightlines. We will discuss the 
results presented here in Section r5.2. 21 

4.3. A Supernova- Driven Interstellar Medium 

In this section, we examine another model in which the 
interstellar gas is heated by SN explosions. This model 
is distinct from the previous model in a number of ways. 
The previous model considered only isolated SNRs above 
z = 130 pc, and the X-ray spectra were calculated from 
ID hydrodynamical simulations of individual remnants. 
Here, we use a 3D hyrodynamical simulation of verti- 
cally stratified interstellar gas that is heated by discrete 
SN explosions (Joung &: Mac Low.,2Q06.) . Unlike the pre- 
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vious model, some SNe occur in clusters, the evolving 
SNRs can interact, older remnants may be re-energized 
by new SNe, and SNe in the Gal actic disk drive a. foun- 
tain of hot gas up into the halo ([Shapiro fc Field! [l976l : 
lBregmanl ll980\ However, ionization equilibrium was as- 
sumed, and magnetic fields neglected. 

4.3.1. Model Description 

The hydrodynamical siniulation used here is described 
fully in iJoung fc Mac Lowl ()2006[ ). and the reader is re- 
ferred to that paper for more details. The simulation 
was carried out using FlaslFl. a parallelized Eulerian hy- 
drodynamical code with adaptive mesh refinement. The 
simulation box extends from z — —5 kpc to z = -1-5 kpc, 
and has a size of (0.5 kpc)^ in the xy plane. The maxi- 
mum spatial resolution is 1.95 pc. The upper and lower 
boundaries have outflow boundary conditions, while the 
vertical sides of the simulation box have periodic bound- 
ary conditions. 

The simulation box was initialized with 1.1 x 10** K gas 
in hydrostatic equilibrium. This gas was then heated by 
discrete SN explosions, each of which injected 10^^ erg of 
energy into a small region of the grid. These explosions 
generally occurred randomly in time and space, with a 
rate 

r{z)^ria,{z) -l-rii(z) 

= (6.2e-l^l/325 P<= + 167e-l^l/^° p'^) kpc"^ Myr-\ (7) 

although 3/5 of the Type H SNe occurred in clusters of 
7 to ~40 explosions. The above rat es differ from those 
in Equ ation ^ in two ways. Firstly, IJoung fc Mac Low 
( 20061) assumed higher Galactic SN rates than iFerriere 
(fT998l) : 1/330 yr"! versus 1/445 yr"! (Type I) and 
1/44 yr^^ versus 1/52 yr~^ (Type II). Secondly, Equa- 
tion (J6l) considers only isolated Type II SNe, whose aver- 
age height is larger than the a verage height o f all Type II 
progenitors (266 versus 90 pc: lFerrierelll995[ ). The model 
also includes diffuse heating, due to photoelectric emis- 
sion from UV-irradiated dust grains, and radiative cool- 
ing. Figure [9] shows vertical slices of the density and 
temperature a.t t = 120.0 Myr. 

We calculated X-ray spectra for 242 sightlines look- 
ing vertically upward and downward from the Galactic 
midplane. The viewpoints of these sightlines formed an 
11x11 grid in the xy plane, with a grid spacing of »49 pc. 
The sp ectra were calculated using the Ravmond fc SmithI 
(|1977l ) spectral code (updated by J. C. Raymond & B. W. 
Smith, 1993, private communication with R. J. Edgar), 
assuming that the material along the line of sight is in 
CIE. We assumed that the gas on the grid is optically 
thin. We ignored the emission from the first 100 pc 
of each sightline, as such material is not in the Galac- 
tic halo. In our observational analysis, emission from 
within ^100 pc of the midplane is attributed to our fore- 
ground model component, derived from ROSAT shad- 
owing data. 

As with the previous model, we used the method de- 
scribed in Section [4.2.21 to characterize the model spec- 
tra with IT models, and used the resulting X-ray tem- 
peratures and X-ray emission measures to compare with 

^^ Developed at the University of Chicago Center for Astrophys- 
ical Thermonuclear Flashes; http://flash.uchicago.edu/web/ 
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Figure 9. Vertical slices of the density (left) and temperature 
(right) at t = 120.0 M yr, from the SN-driven ISM model of 



IJoung fc Mac Low! 1 120061) . Note that the full simulation box ex- 
tends to l^l = 5 kpc. 



our observations. As in Section I4.2.2[ our simulated 
XMM-Newton spectra included foreground, extragalac- 
tic background, and instrumental background compo- 
nents, as well as the halo emission from the SN-driven 
ISM model. The emission from the halo and extra- 
galactic components was subjected to absorption with 
TVh — 1-7 X 10^*^ cm~^. This column density was also 
used in the subsequent spectral fitting, from which we ob- 
tained the characteristic X-ray temperatures and X-ray 
emission measures. Figure [TU] shows a simulated XMM- 
Newton spectrum for one of our model sightlines, along 
with the best-fitting model. The differential emission 
measure for the input halo spectrum is shown in Fig- 
ure [7| (dashed line). Because we assume that the plasma 
is in CIE, only plasma with T > 10^ K will contribute 
in the XMM-Newton band. Although the input model 
predicts emission from a range of temperatures, a model 
with a IT halo generally fits the simulated spectra well. 
It should be noted that our model spectra do not 
take into account absorption by cold gas on the grid - 
when we characterized our model spectra, we assumed 
that the absorbing column (N^ = 1.7 x 10^° cm~^) 
was located entirely beneath the hot X-ray-emitting gas. 
For ^2/3 of the model sightlines, the column density 
of cold gas that is mixed in with the bulk of the X- 
ray-emitting gas is <10^^ cm~^ (i.e., an order of mag- 
nitude smaller than the column densities used in our 
XMM-Newton analysis). We have investigated the ef- 
fect of ignoring on-grid absorption by creating simulated 
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Figure 10. As Figure [6] but using the SN-driven ISM model as 
the input halo model. 

XMM-Newton spectra from an XSPEC model of the 
form phabs * (raymond + phabs * raymond), and char- 
acterizing the resulting spectra with a model of the form 
phabs * raymond. As with our previous simulations, the 
column density of the first phabs component was fixed 
at 1.7 X 10^° cm~^. The two raymond components in the 
input model had the same emission measure, and temper- 
atures of 1.5 X 10^ K and 2.5 x 10^ K (it does not matter 
which component is the hotter - our conclusion is the 
same cither way). We found that the characteristic X- 
ray temperatures and X-ray emission measures obtained 
were not strongly affected by the column density between 
the two raymond components in the input model, at least 
up to column densities of a few x 10^" cm~^. We there- 
fore conclude that ignoring on-grid absorption will not 
adversely affect our results. 

In the following, we also compare the predicted X-ray 
surface brightnesses with our observations. These values 
were extracted directly from the model spectra. 

4.3.2. Comparing the SN-Driven ISM Model with 
Observations 

Figure [11] shows the time variation of the X-ray spec- 
tral properties predicted by the SN-driven ISM model. 
The model has not settled down to a steady state - the 
typical X-ray temperature and X-ray surface brightness 
rise steadily throughout the period shown. However, in 
this temperature regime, an increase in X-ray surface 
brightness can be brought about by an increase in tem- 
perature as well as by an increase in emission measure. 
The typical X-ray emission measure does not rise steadily 
throughout the plotted period. 

Figure [TT] also compares the predicted X-ray proper- 
ties with the observed values from our XMM-Newton 
analysis. Similarly to Section 14.2.31 we multiplied our 
observed emission measures and surface brightnesses by 
sin \b\ before comparing them with the model values. At 
the earlier times, the model halo is too cool and too faint 
in the 0.4-2.0 keV band. Around t = 110 Myr, the pre- 
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Figure 11. Comparison of the halo (a) X-ray temeratures, (b) X- 
ray emission measures, and (c) intrinsic 0.4-2.0 keV surface bright- 
ness predicted by the SN-driven ISM model with the values ob- 
tained from our XMM-Newton observations. The boxplots show 
the predicted values for a range of simulation times. The rightmost 
datapoint shows the results obtained by averaging the spectra from 
several of the timesteps (see text for details). Each box indicates 
the median and quartiles, and the whiskers extend to the most ex- 
treme point that is no more than 5 times the interquartile range 
from the box; any outliers are plotted individually with open cir- 
cles. The horizontal dashed lines indicate the observed median 
and quartiles, and the gray band indicates the 90% bootstrap con- 
fidence interval on the observed median. 

dieted X-ray temperatures are in reasonable agreement 
with the observed value, but the halo is typically a factor 
of ^4 too faint. At later times the model halo is too hot; 
by t = 155 Myr it is also somewhat brighter than is ob- 
served, although the median predicted surface brightness 
is within 50% of the observed median. 

From Figure [TT] is it not clear what state the h alo is 
tending toward HOe Avillez fc BreitschwerdO ([2004[ ) esti- 
mated that a steady state halo in a simulation such as 
this should be reached in < 180 Myr. However, it is possi- 
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ble that, instead of settling down to a steady state that is 
hotter and brighter than the observed halo, the variation 
in X-ray temperature and surface brightness is part of a 
slow oscillation about a mean state, with a period of at 
least ^^130 Myr. Determining the final state predicted by 
this model would probably require running it for at least 
several more tens of megayears, which is not currently 
practical. We therefore try two different approaches in 
comparing the simulation data to our observations. 

Our first approach is to assume that the data shown 
in Figure [11] represent roughly half a cycle in the oscil- 
lation of the halo about a mean state. The simulation 
models only a narrow column of the halo. As our observa- 
tion directions are not all toward the Galactic poles, our 
sightlines sample different spatial locations in the halo, 
which, in a statistical sense, should correspond to differ- 
ent times throughout the cycle. We therefore proceed by 
averaging the spectra for each sightline from 7 of the 8 
timesteps shown in Figure[TT](we don't use the data from 
i = 118 Myr, so as not to oversample the times around 
t ^ 120 Myr). We characterize the resulting averaged 
spectra with IT models, using the procedure described 
previously. The results are compared with our obser- 
vations in Figure [11] (the rightmost datapoint, labeled 
"Ave"), and in Figures [T^ andlT^. 

Our second approach is to assume that the final time 
in Figure [TT] ft = 155 Myr) gives our best estimate of the 
steady state that the halo is approaching. We use the 
predictions from that final time in Figures [T2h and[T2H. 

Whether we use the averaged spectra or the final 
timestep, the predicted halo temperature is too high. 
The median predicted temperatures are 2.73 x 10^ K 
(averaged spectra) and 3.55 x 10^ K (final timestep), 
against an observed median of (2.08^0 07) ^ 10^ K. The 
averaged spectra underpredict the halo emission mea- 
sure, although the discrepancy is not as large as with 
the extraplanar SNR model. The median predicted X- 
ray emission measure is 0.44^q'52 dex smaller than the 
median observed value of E.M. sin |&|: 0.69 x 10"'^ versus 
(1.90^046) ^ 10~^ cm~^ pc. For comparison, the dis- 
crepancy for the extraplanar SNR model is 1.371;q'52 ^ex. 
For the spectra from the final timestep, the median pre- 
dicted X-ray emission measure (1.17 x 10~^ cm~^ pc) is 
in better agreement with the observed median (the dis- 
crepancy is 0.2lj:°;52 dex). The predicted 0.4-2.0 keV 
surface brightnesses are generally in reasonably good 
agreement with the observed values, although there are 
a few sightlines in the final timestep of the model that 
exhibit much greater surface brightnesses. The me- 
dian predicted values are within 50% of the median 
observed value of 5'o. 4-2.0 sin |fe|: 1.04 x 10^^^ (aver- 
aged spectra) and 2.37 x 10^^^ (final timestep), versus 
(1.62to.o5) X 10^^^ erg cm^^ g^i deg"^ (observ ed) W e 
will discuss the results presented here in Section [5.2.31 

5. DISCUSSION 

5.1. The Temperature and Emission Measure of the 
Galactic Halo 



other recent studies. 

Our measured halo temperatures are typically ^2 x 
10^ K - the median temperature is 2.08 x 10® K, and the 
lower and upper quartiles are 1.96 x 10® and 2.24 x 10® K, 
respectively. The halo temperatures measured here, us- 
ing a thermal plasma model, are more tightly constrained 
than those inferred from the O vil/O Vili intensity ratio 
(Paper I). This is partly because fitting with a thermal 
plasma model uses more of the information in an ob- 
served spectrum, and partly because the errors on the 
intensities of individual lines are not combined in the 
temperature uncertainty. 

Figure 113b . compares our halo temperatures with 
various other measurement s made with ROSAT , 
XMM-Newton^ and Suzak u (iKuntz fc S nowden' 120001: 
Smith et all [20 07: Galc azzi et all 120071: Yoshino et all 
120091: iLei etaJI 120091: IGupta et all I2009D . For papers 
with multiple results, the figure caption indicates the 
specific results that we have plotted. Our temperatures 
are typically lower than the temperature obtained by 
IKuntz fc Sno^dCTJ ([20Q0l ). using ROSAT All-Sky Sur- 
vey data. The most likely reason for this is that 
IKuntz fc Snowdeiil (2000) used a two-temperature (2T) 
model for the halo, whereas as we used a IT model. An 
additional, lower-temperature halo component could ac- 
count for some of the lower-energy flux in the XMM- 
Newton band, enabling the other component to shift 
toward higher temperatures. However, in practice we 
found that we were unable to constrain a 2T model us- 
ing our XMM-Newton data, as XMM-Newton^s sensitiv- 
ity does not extend to as low energies as RO SAT's. A 
IT halo model adequately described our XMM-Newton 
spectra. 

Our halo temperatures are gen erally in goo d agree- 
ment with those measured bv iSmithetall (I2007D. 
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iGaleazzi et"aI1 ([200l. 'Gupta etal.1 (l2009ri. 
(|2009D . and [Yoshino et'al.. (,2009[ l All of these mea- 
surements are from IT halo models. Note that the 
I Yoshino et al.l ([2009 !) temperatures were obtained with 
a fixed foreground model for all their sightlines, and Fe 
and Ne abundances (relative to O) that were free to 
vary (their Table 6). When they fitted their model with 
an independent f oreground mod el and fixed abundances 
for each sightline. [Yoshino et al.l obtained systematically 
high er halo temper atures (their Tables 3 and 4). How- 
ever, [^shmoler^l] argue that these higher temperatures 
are mostly determined by the Fe-L and Ne-K emission, 
and are inconsistent with the behavior of the O vii and 
O VIII emission. When they measured th e halo temper- 
atures that we have used in Figure [TSk . [Yoshino et al.l 
found that about half of their sightlines required either 
an overabundance of Fe and Ne, or an additional hotter 
emission component. 

The halo temperatures measured here and in other re- 
cent studies are fairly constant across the sky. This is 
most likely because gas with T < 1 x 10® K would be dif- 
ficult to detect with XMM-Newton (unless it has a large 
emission measure), while gas with T > 3 x 10® K would 
escap e from the Galactic potential well (e.g., [Bregmaiil 



In Section [521 below, we will discuss the various phys- 
ical models of the hot halo that we examined in the pre- 
vious section. Here, we discuss our XMM-Newton halo 
measurements, and compare the results with those from 



Figure [T3b compares our halo emission measures with 
those from the studies discussed above. In general, our 
values are in good agreement with the values from the 
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Figure 12. Histograms comparing tiie iialo X-ray temperatures (panel a and c) and X-ray emission measures (panels b and d) predicted 
by the SN-driven ISM model (solid gray) with those obtai ned from our XMM-Newton observations (black outline). The observed emission 
measures have been multiplied by sin |6| (see Secti on | 4.2.3l l. The model temperatures and emission measures were derived by characterizing 
the model spectra with IT models (see Section I4.2.2| l. The model histograms have been rescaled so they cover the same area as the 
correspo ndin g observational histograms. Panels (a) and (b) show the results obtained by averaging the spectra from several model tim esteps 
in Figure [TT] (see text for details), while panels (c) and (d) show the predictions from the final timestep (i = 155 Myr) in Figure [TT] 



oth er studies . Th e large emission measure obtained 
by 'Lei et al.' ('2009') may be partly due to their us- 
ing Wilms ct al. (2000) abundances. The I Wilms et al.l 
(120001) oxygen abundance is a factor of 1.8 smaller than 
thel Anders fc Grevessd (|1989D value. As oxygen emission 
tends to dominate the halo X-ray spectrum, a smaller 
oxygen abundance will result in a larger X-ray emission 
measure. 

There is considerable scatter in the halo emission mea- 
sure, with the values spanning an order of magnitude 
(~0.0005-0.006 cm-6 pc). This patchiness to the halo 
has already been pointed out by lYoshino et al.l (|2OO90 
and in Paper I. We will discuss in more detail the various 
physical models that we have examined in Section 15.21 
below, but we note here that a patchy halo favors an 
inhomogeneous, stochastic heat source, such as SNe, as 
opposed to accretion of extragalactic material, which we 
would expect to be fairly homogeneous. 

5.2. Physical Models of the Hot Halo 

Here, we discuss each of the models presented in Sec- 
tion |4] in turn. 

5.2.1. The Disk Galaxy Formation Model 

In Se ction |4. II we examined a disk galaxy fo rmation 
model (|Toft et al.l [20021: iRasmussen et al.ll2009[ ). which 
predicts the existence of a hot halo extended over tens of 
kiloparsecs, formed from material falling into the galaxy's 



poten tial well. Using results from IRasmussen et al.] 
(|2009D . we find that the emission- weighted mean tem- 
perature predicted for a Milky Way-sized galaxy (vc ^ 
220 km s~^) is in good agreement with the median ob- 
served halo temperature, but the model underpredicts 
the X-ray luminosity of the halo by at least an order of 
magnitude. 

The above results suggest that the extended hot halo 
predicted by disk galaxy formation models is not a major 
contributor to the halo X-ray emission that we observe 
with XMM-Newton. However, this interpre t ation is com- 
plicated by the simulations of iCrain et al.l (I2010D. Thev 
argue that the stellar feedback in the IRasmussen et al.l 
()2009() simulations is too weak, resulting in too much 
mass ending up in stars and too little mass in hot gas. 
Ideally, we would compare the surfa ce brightness of the 
accreted extragalactic material in the ICrain et al.l (|2010[ ) 
simulations with the halo surface brightness obtained 
from our XMM-Newton observations, but this is not pos- 
sible for two reasons. Firstly, in addition to inf ailing 
extragalactic material, the ICrain et al.l model includes 
interstellar gas that was heated by stars and SNe and 
transferred upwards by an outflow (similar to the galac- 
tic fountain in the model described in Section H75t . Sec- 
ondly, their model predictions are in terms of luminosity 
rather than surface brightness. Because gas heated by 
stars and SNe will tend to be concentrated nearer the 
disk, compared with the more extended halo of accreted 
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Figure 13. Halo (a) temperatures and (b) emission measures 
plotted against Galactic latitude. The open circles show the results 
from this work, with the erro r bars omitted for cla rity. The solid 
triangles show the results from i Yoshino et al.l 1)20091 ) with |6| > 20° ; 
specifically, the results obtained with a fixed foreground model for 
all sightlines (their Table 6). The solid square shows t he results 
from a Suzaku shadowing observation of MBM 12 (S mith et al.l 
120071 ) . The "+" and the " X " show the results from shadowing ob- 
servations of MBM 20, carried out with XMM-Newton and Suzaku, 
respectively (Galeazzi et al. 2007; Gupta et al. 2009). Specifically, 
the XMM-Newton result is the result obtained with frozen abun- 
dances from Table 2 in Galeazzi ct al. (2007), and the Suzaku result 
is from row 3 of Table 2 in Gupt a et al. ( 2009|) . The star shows the 
results from a Suzaku shadowing observa tion of a d usty filament 
in the southern Galactic hemisphere (Lei ct al."2009'); specifically, 
model 2 from their Table 1. The horizontal dashed line shows the 
temperature derived bv lKuntz fc Snowden. (.2000, ) from the ROSAT 
All-Sky Survey. 

extragalactic material, it may provide a much larger sur- 
face brightness per unit luminosity than the extended 
halo gas. Predictions of the separate contributions to 
the surface brightness due to the gas heated by stellar 
processes and the extended halo are not currently avail- 
able from the iCrain et al.l model. Such predictions are 
needed to determine whether or not accreted extragalac- 
tic material makes a major contribution to the observed 
halo X-ray emission. 

Although surface brightness predictions are not cur- 
rently available, we can make rough comparis ons between 
the t otal X-ray luminosity predicted by ICrain et al.l 
(|2010f ) and the halo luminosity expected from the XMM- 
Newton data. An Lx-T plot derived from their simu- 
lations shows that galaxies with emission-weighted halo 
temperatures of ^2 x 10^ K have 0.3-2.0 keV luminosities 
of ~3-40 X 1039 erg s^^ (R. Grain, 2010, private com- 
munication) , compared with <0.2 x 10^^ erg s^^ from 
IRasmussen et al.l (|2009[ ). In Section 14.11 we derived a 
Galactic halo luminosity of 1.9 x 10^^ erg s~^ from our 



observations, which is sli ghtly smaller than the range of 
luminosities predicted bv iGrain et al.l (|2010D . However, 
this observed luminosity was calculated assuming that 
the observed emission came from a uniform sp h ere o f 
radius i?sph = 15 kpc. Figure 2 in IGrain et al.l ()2010f ) 
implies that the emission in their model is much more 
extended than this, with i?sph ~ 50 kpc. Using this 
larger radius leads to an observed halo luminosity of 
1.3 x 10^° erg s ~^, which lies within the range predicted 
bv IGrain et al.l . 

5.2.2. The Extraplanar SNR Model 



In Section 14.2! we examined a model in which the hot 
halo gas is co ntained i n an ensemble of isolated extra- 
planar SNRs (|Shelton|[2006[) . The X-ray temperatures 
predicted by this model are in good agreement with the 
observed values. However, the predicted X-ray emission 
measures, obtained from model spectra above 0.4 keV, 
are typically an order of magnitude smaller than the ob- 
served values. We can increase the predicted X-ray emis- 
sion measures by increasing the SN explosion energy, Eq , 
the effective ambient magnetic field, B^s, or the SN rate. 
However, doubling all three of these parameters still re- 
sulted in X-ray emission measures that were too small. 

Of the above three parameters, B^s has the largest 
effect on the predicted X-ray emission measures (see Ta- 
ble HJ . Increasing B^s increases the non-thermal pres- 
sure, and therefore increases the compression of the hot 
X-ray-producing bubble. This increased compression in- 
creases the density and temperature of the bubble, in- 
creasing the X-ray emission measure inferred from the 
predicted emission in the XMM-Newton band. We have 
not carried out simulations with Bcff > 5.0 fiG, but if 
we extrapolate the values in Table |4l we find we would 
need a B^s of a few tens of microgauss to match the ob- 
served emission measures. Such a magnetic field implies 
a non-thermal pressure Pnt/k > 10^ cm"^ K. This is 
an implausibly high non-thermal pressure for the halo, 
as it is several times larger than th e midplane value 
(jBoulares fc Goxlll990t [Ferrierell200lD . In addition, the 
model with Bcs = 2.5 /iG predicts that hot gas would be 
seen on only ~l/3 of sightlines. As noted above, increas- 
ing Beg results in smaller, hotter remnants. In addition, 
these smaller remnants are brighter and so shorter-lived. 
As a result, increasing BcS means that even fewer sight- 
lines would intercept remnants. In reality, we find hot 
gas on most, if not all, of our XMM-Newton sightlines. 
Therefore, we cannot bring the extraplanar SNR model 
into agreement with our observations by increasing the 
assumed effective ambient magnetic field. 

Increasing the SN rate increases the number of SNRs 
expected to lie along a given sightline, and so increases 
the predicted X-ray emission measures. However, we 
would have to increase the SN rate given by Equation ([6]) 
by a factor of ~6 in order to give the same increase in 
the median X-ray emission measure that we see when we 
increase B^s from 2.5 to 5.0 /zG. To match the observa- 
tions, we would have to increase the SN rate by an even 
larger factor. The Galactic SN rate of ^2 per century is 
constrained to within a factor of ~2 (see the online Sup- 
plementary Informatiorn for iDiehl et al.l I2006D . If the 
SN scale heights in Equation (j6|) are well constrained, 

^^ http://www.nature.com/nature/journal/v439/n7072/suppinfo/nature0436' 
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the halo SN rate is also constrained to within a factor 
of ~2. Even if we double the SN scale heights in Equa- 
tion dni, the integrated SN rate above \z\ = 130 pc only 
increases by a factor of 2.5. Therefore, to bring the extra- 
planar SNR model into agreement with our observations 
we would have to increase the SN rate by an unrealistic 
amount. We also cannot bring the model into agreement 
with our observations by increasing i?o, as Table |4] shows 
that the median X-ray emission measure depends only 
weakly on Eq. 

We therefore conclude that the hot halo gas that we 
observe with XMM-Newton is not primarily due to an en- 
semble of isolated extraplanar SNRs. Most of this popu- 
lation of remnants would be relatively old (age > 1 Myr) 
and faint in the XMM-Newton band. It should be noted, 
however, that this result does not imply that extraplanar 
SNRs do not contribute to the hot halo gas at all. In the 
XMM-Newton band (above 0.4 keV), their contribution 
is masked by brighter emission from an SN-driven galac- 
tic fountain (see Section [5.2.3[ below). At lower energies, 
extraplanar SNRs could still contribute significantly to 
the hot gas observed in the 1/4 keV band, as originally 
suggested by Shelton (2006). In addition, young, bright 
remnants, although r are, do produce eniission that is de- 
tectable by Suzaku (jHenlev fc Sheltonl I2OO90 and that 
should also be detectable by XMM-Newton. 

5.2.3. The SN-Driven ISM Model 



3.5 



In Section 14.31 we examined a hydrodynamical simu- 
lation of vertic ally stratified interstella. r gas, driven by 
SN explosions (iJoung fc Mac Lowll2006l ). Unlike the ex- 
traplanar SNR model, the SNRs do not evolve in iso- 
lation, and the model r esults in a galactic fountain of 
hot gas up into the halo ([Shapiro fc Fieldlll976l : iBregmaiil 
[1980). Also, for this model we assumed that the gas is in 
CIE, whereas the extraplanar SNR model includes self- 
consistent modeling of the ionization evolution. 

As with the extraplanar SNR model, we folded the 
model spectra through the XMM-Newton response, 
added photon noise, and characterized the resulting spec- 
tra with IT models. Before going on to discuss the com- 
parison of the resulting X-ray temperatures and emission 
measures with our XMM-Newton observations, we shall 
look at how well these X-ray spectral properties reflect 
the properties of the gas on the hydrodynamical grid. 

Figure 114b . compares the X-ray temperatures derived 
from the IT models with the mean temperature along 
each model sightline through the hydrodynamical grid, 
weighted by the 0.4-2.0 keV emission. We find that the 
X-ray temperatures are in reasonable agreement with the 
emission-weighted mean temperatures (typically within 
0.4 X 10^ K). Figure fT4b compares the X-ray emission 
measure for each sightline with the emission measure of 
the gas with T > 10^ K along that sightline (we use 
this quantity because there is no obvious analog to the 
emission- weighted mean temperature). The X-ray emis- 
sion measure underestimates the emission measure of gas 
with T > 10^ K, although the two quantities are well cor- 
related and generally agree within a factor of 2.5. Over- 
all, we conclude that the properties derived from the IT 
models are good characterizations of the hot gas on the 
hydrodynamical grid 

The 0.4-2.0 keV halo surface brightness predicted by 
this model is in good agreement with the observations 
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Figure 14. Comparison of X-ray spectral properties derived from 
IT models (ordinates) and properties derived directly from the 
hydrodynamical data (abscissae) for each sightline through the SN- 
driven ISM model at t = 120 Myr. 

(within ^50%). As the emission predicted by this model 
is much brighter, it will mask the contribution from ex- 
traplanar SNRs. Note that this model does include some 
in situ heating by extraplanar SNe. The integrated halo 
SN rate above \z\ ~ 100 pc is ~2 times larger in this 
model than in the extraplanar SNR model (although the 
volumetric SN rate falls off more rapidly with |z| in this 
model). However, this larger integrated SN rate is not 
enough to explain the differences in X-ray surface bright- 
ness between the models: this model predicts emission an 
order of magnitude brighter than the extraplanar SNR 
model. The majority of the halo X-ray emission in this 
model comes from hot gas that is driven from the disk 
into the halo by a galactic fountain. As mentioned in Sec- 
tion [5T231 it is possible that the emission from a galactic 
outflow or fountain also masks the contribution from an 
extended hot halo of accreted material. 

Although the median predicted X-ray surface bright- 
ness and X-ray emission measure are in good agreement 
with the observed medians, the halo X-ray temperature 
predicted by this model is ~1 x 10^ K larger than our 
measured value, independent of whether the simulated 
halo is undergoing a long-period (P > 100 Myr) oscil- 
lation about some mean state, or is settling down to a 
steady state. Given the uncertainties in the model pre- 
dictions (due to the time variability shown in Figure [TT]), 
these discrepancies are probably insufficient to rule out 
this model. Nevertheless, we discuss below three possible 
causes for the temperature discrepancy: (1) the simula- 
tion overpredicts the temperature of the hot gas in the 
halo; (2) the temperature of the hot gas does not ac- 
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curately predict the X-ray spectrum, because the gas is 
not in ionization equilibrium; or (3) a bias in our spec- 
tral analysis causes us to underestimate the observed halo 
X-ray temperature. Correcting for whichever of these ef- 
fects turn out to be important may also affect the X-ray 
surface brightnesses and emission measures. However, 
it is not easy to foresee whether these corrections will 
increase or decrease the discrepancies between the pre- 
dicted and observed surface brightnesses and emission 
measures. 

The first possible cause of the temperature discrep- 
ancy is that the simulation overpredicts the temperature 
of the hot gas in the halo. The simulation does not in- 
clude thermal conduction, which could in principle lower 
the temperature of the hot X-ray-emitting gas, although 
thermal conduct ion in the ISM can be sup p ressed by 
magnetic fields. 'De Avillcz fc Brcitschw erdtl ()2004| ) ar- 
gue that turbulent diffusion is more cfhcient at mixing 
hot and cold gas than thermal conduction. However, the 
efficiency of this mixing will be underestimated if the 
mixing is not fully resolved. In the simulation used here, 
the spatial resolution is between 1.95 and 15.6 pc, but 
the adaptive mesh refinement criterion is chosen to focus 
maximum numerical resolution on the 200 pc above and 
below the midplane, so the gas is less well resolved in 
the halo than in the disk. As a result, it is possible that 
turbulent mixing of hot and cold gas is under-resolved in 
the halo, resulting in an ov erestimate of the t emperature 
of the X-ray-emitting gas (Fuiit a et al.ll2009r ). 

A lack of spatial resolution may also suppress the radia- 
tive cooling of the hot gas - averaging the hot gas density 
over large cells eliminates local denser reg ions that would 
radia t e more efficiently (de Avillcz fc Breitschwerdtl 
120041 ). iDe Avillez k Breitschwerdtl (|2004i) found that the 
filling factor of hot (T > 10^'^ K) gas was not signif- 
icantly affected by the spatial resolution of their sim- 
ulations. However, there is insufficient information to 
determine how the emission-weighted mean temperature 
of the X-ray-emitting gas would be affected. Simulations 
with higher spatial resolution in the halo will help deter- 
mine whether or not the mixing of hot and cold gas and 
the radiative cooling of hot gas are adequately resolved 
fo r predicting the X- r ay em ission. 

iJoung fc Mac Low! (|2006l ) pointed out that the aver- 
age gas density at several disk scale heights and beyond 
(0.2 kpc < \z\ < 2.5 kpc) is somewhat underpredicted in 
their model compared to observations. They suggested 
that additional components of pressure from the mag- 
netic field and cosmic rays may contribute significantly 
to the support in the vertical direction (see Section 3.1 
of their paper). Additional vertical pressure support will 
lead to larger disk scale heights and hence larger gas 
masses and lower temperatures at 1-2 kpc heights, which 
will reduce the discrepancy between the observations and 
the SN-driven ISM model. Preliminary results from mag- 
netohydrodynamics simulations are in agreement with 
this expectation (A. Hill et al., 2010, in preparation). 

The second possible cause for the discrepancy between 
the predicted and observed X-ray temperatures is that 
the hot gas is out of equilibrium. This gas was shock- 
heated by SNe - this rapid heating causes the ionization 
temperature (which determines the X-ray spectrum) to 
lag behind the kinetic temperature (which is the quan- 



tity obtained from the hydrodynamical simulation). This 
ionization evolution was followed self-consistently in the 
ID SNR models described in Section 14. 2[ but not in 
this 3D simulation. CIE is reached on a timescale of 
tpn ^ lO^^n^ ^ s, where rie is the electron density in cm~'^ 
(jMasailll994[ ). In the simulation, the density of gas with 
T > 10** K is <10-3 cm-3, implying igq > 30 Myr. This 
is similar to the dynamical timescale (5 kpc/cs ~ 30 Myr, 
where Cg is the sound speed in 10® K gas), and so we 
would expect the hot gas to be at least partially underi- 
onized. 

Calculating the X-ray spectrum of an underionized 
plasma requires knowledge of the ionization balance, 
which in turn depends on the history of the plasma. How- 
ever, we note that an underionized plasma will have fewer 
higher-stage ions (e.g., 0+^) relative to lower-stage ions 
(e.g., O"*"®) than we would expect from the kinetic tem- 
perature. Therefore, if the halo is underionized, mod- 
eling the observed emission with a CIE model will un- 
derestimate the kinetic temperature of the halo gas. We 
have confirmed this by simulating XMM-Newton obser- 
vations of an underionized plasma (using the XSPEC nei 
model) , and characterizing the resulting simulated spec- 
tra with IT CIE models (for simplicity, here we just use 
the nei model as an input for our simulations, instead 
of the full multicomponent model used in Section 14.2.2^ . 
For Tipt < 10^ cm~'^ s, CIE models do not fit the simu- 
lated spectra well (these models underestimate the low- 
energy flux). However, for rict ^ 10^-10^^ cm~'^ s, IT 
CIE models give good fits to the simulated spectra, but 
underestimate the input kinetic temperature. 

The above discussion considers only the hot gas be- 
ing underionized. An additional possible source of X- 
rays is delayed recombination from overionized gas - gas 
that has undergone rap id adiabatic co oling as it expands 
into the halo ( Breitsch werdt fc Schm utzler 1994, 1999). 
With our current simulation data, it is difficult to esti- 
mate the contribution of delayed recombination to the 
halo emission, relative to the emission from the hot gas. 
Simulations that can track the ionization evolution of 
the plasma are needed to determine the extent to which 
the assumption of CIE affects the predicted spectra, and 
hence the derived X-ray temperature. 

The third possible cause of the temperature discrep- 
ancy is that something (other than the assumption of 
CIE) is biasing the temperatures measured from our 
XMM-Newton observations. The discrepancy is not be- 
cause we characterize the halo e mission with a IT model , 
as opposed to a 2T model (e.g.. lKuntz fc S nowden 20001: 
see Section IO]) . because we characterize the model spec- 
tra in the same way. If we are underestimating the halo 
temperature, the most likely cause is that our foreground 
(LB and/or SWCX) model is too faint; in particular, that 
it underestimates the foreground O vii emission (under- 
estimating the foreground O vii emission means that the 
halo O VII emission is overestimated relative to the halo 
O VIII emission, causing the measured halo temperature 
to shift to a lower value). However, our halo tempera- 
tures are in good agreement with those from other recent 
studies (see Section ISTTj) . which used different methods to 
estimate the foreground emission. It therefore seems un- 
likely that a bias in the halo temperature measurements 
is causing the discrepancy between the observed temper- 
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atures and those predicted by the SN-driven ISM modeL 
Nevertheless, an accurate model of SWCX emission will 
help evaluate our method for estimating the foreground 
emission. 

6. SUMMARY AND CONCLUSIONS 

We have analyzed 26 high-latitude XMM-Newton ob- 
servations of the SXRB, concentrating in particular on 
the emission from the Galactic halo. These observations 
were chosen from a much larger set of observations (Pa- 
per I) as they are expected to be the least contaminated 
by SWCX emission. We modeled the 0.4-5.0 keV X-ray 
spectra with emission components from the foreground, 
the Galactic halo, and the extragalactic background, 
with additional components modeling parts of the in- 
strumental background. Assuming a single-temperature 
CIE plasma model for the halo, we typically obtained 
halo temperatures of ~(1.8-2.4) x 10^ K, and emission 
measures of ^0.0005-0.006 cm~^ pc, in good agreement 
with previous studies. While the halo temperature is 
fairly constant, the emission measure exhibits significant 
sightline-to-sightline variation. 

We compared the observed X-ray properties of the 
halo with the predictions of three physical models for 
the origin of the hot halo gas: (1) a disk galaxy for- 
mation model, which predicts the ex istence of a hot 
halo extended over tens of kilop arsccs ("T oft et all 120021 : 
[Rasmusscn ct al. 2009; Grain et al. 2010); (2) a model in 
which the halo is heated by extraplanar S Ne, an d the 
hot gas resides in isolated SNRs (SheltonI 120061 ) ; and 
(3) a more comprehensive model of SN heating of the 
ISM, in which the SNRs do not evolve in isolation, but 
drive a fountain of hot g as from the disk into the halo 
(jJoung fc Mac Lowll200l . 

Model 2 matches the observed halo temperature rea- 
sonably well, but this model predicts emission at least an 
order of magnitude too faint in the XMM-Newton band, 
implying that another source of hot gas is needed, in 
addition to isolated extraplanar SNRs. With Model 1, 
the conclusions are more uncertain: the original simu- 



lations that we examined (from iRasmussen et al.1 1200£ 
predicted luminosities at least an order of magnitude too 
faint in the XMM-Newton band, whereas newer simula- 
tions with st ronger stellar feed back predict larger X-ray 
luminosities (jCrain et al.ll20ldl ) , in better agreement with 
those inferred from our observations. However, the pre- 
dicted emission from that model includes contributions 
not only from accreted extragalactic material falling into 
the galactic potential well, but also from material that 
has been heated by stars and flowed out from the disk. 
Predictions of the relative X-ray surface brightnesses due 
to these two processes are needed to determine whether 
or not emission from accreted extragalactic material is a 
major contributor to the observed halo X-ray emission. 
A flow of hot gas from the disk into the halo, in the 
form of an SN-driven galactic fountain (jShapiro fc Fieldl 
119761: iBregmanI [19801) ■ also occurs in the SN-driven ISM 
model (Model 3). The X-ray surface brightness predicted 
by this model is in good agreement with the observed 
surface brightness of the halo, although it should be 
noted that the halo in this model has not yet settled 
down to a steady state. Therefore, while we cannot cur- 
rently rule out the possibility of a significant contribution 
from accreted extragalactic material, our analysis indi- 



cates that hot gas in an SN-driven galactic fountain is a 
major (possibly dominant) contributor to the halo X-ray 
emission. While previous work has shown that di sk SNe 



produ ce sufficient energy to heat the ha lo Ce.g.. iWane 
1991 iShelton et aD[2007t lYao et al.l[200l lYoshino et al 



2009( 1. to the best of our knowledge this is the first time 
that CCD-resolution spectra of the halo have been com- 
pared with predictions from a hydrodynamical model of 
a galactic fountain. 

Although the halo X-ray surface brightness predicted 
by the SN-driven ISM model is in good agreement with 
the observations, this model overpredicts the X-ray tem- 
perature of the halo. We discussed various possible 
reasons for this discrepancy, including the simulation 
under-resolving the mixing of hot and cold gas, the sim- 
ulation omitting important sources of pressure support 
that would increase the disk scale height and lower the 
temperature at \z\ ^ 1-2 kpc, the hot halo gas be- 
ing underionized, and our potentially underestimating 
the foreground O VII intensity in the observed spectra, 
which would in turn cause us to underestimate the ob- 
served halo temperature (although this final possibil- 
ity seems unlikely, given the good agreement between 
our temperature measurements and those from other re- 
cent studies). On the modeling side, future simulations 
that have higher resolution in the halo, that include 
magnetic fields, and t hat track the ionization evolution 
(jBenjamin et al.ll2001h may help reduce the discrepancy 
between the predicted and observed X-ray temperature. 
A SWCX model that can more accurately predict the 
foreground emission will help on the observational side. 

We can further test the current and future simulations 
with additional observations. We plan to expand our 
XMM-Newton survey (Paper 1) to cover the whole sky. 
X-ray absorption line measurements can also be used 
to test the models. The significant sightline-to-sightline 
variation in the observed halo emission shows that, in or- 
der to examine halo models in detail, data from as many 
sightlines as possible should be used. 
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